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An evolution of the IR-Radio correlation at very low flux 
densities? 
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ABSTRACT 

In this paper we investigate the radio-MIR correlation at very low flux densities using 
extremely deep 1.4 GHz sub-arcsecond angular resolution MERLIN+VLA observa- 
tions of a 8'.5x8'.5 field centred upon the Hubble Deep Field North, in conjunction 
with Spitzer 24 /im data. From these results the MIR-radio correlation is extended 
to the very faint (~microJy) radio source population. Tentatively we detect a small 
deviation from the correlation at the faintest IR flux densities. We suggest that this 
small observed change in the gradient of the correlation is the result of a suppres- 
sion of the MIR emission in faint star-forming galaxies. This deviation potentially has 
significant implications for using either the MIR or non-thermal radio emission as a 
star- formation tracer of very low luminosity galaxies. 

Key words: galaxies:starburst galaxies: high-redshift infrared: galaxies radio contin- 
uum: galaxies 



1 INTRODUCTION 

Since 1970s and 1980s studies of the radio and far-infrared 
(FIR) properties of galaxies have shown there to be a 
tight correlation between their emission in these two ob- 
serving bands which e xtends over several o rders of mag- 
nitud e in luminosity (Ivan der Kruit 1 1 19731 : ICondon et all 
1982). The advent of the Infrared A stronomical Satel- 
lite (I RAS) All Sky Su rvey in 1983 l|Neugebauer et al\ 
1 19841 : | Soife r et al\ Il987f ) enabled much larger system- 
atic samples of galaxies to be studied at infrared wave- 
lengths and hence further demonstrated the consistency 
and tightness of this corre l ation, albeit for rel a tively 
nearby sources llHelou et all 1 19851: Ide Jong et al\ 1 19851: 



Condon fc Broderickl Il98fj ; ICondon < l ql.lll99ll ; lYun et all 
2001t h Following IRAS, deep observations using the Infrared 



radio correlation holds for relatively bright star-forming 
galaxies (S2pcm > 11 5 /J.3y) out to at least a redshift of 1 
(Appleton et all\2004 ). 

Radio and infrared emission from galaxies in both the 
nearby and distant Universe is thought to arise from pro- 
cesses related to star-formation, hence resulting in the cor- 
relation between these two observing bands. The infrared 
emission is produced from dust heated by photons from 
young stars and the radio emission predominately arises 
from synchrotron radiation produced by the acceleration of 
charged particles from supernovae explosions. It has how- 
ever recently been suggested that at low flux density and 
luminosities there may be some deviation from the tight 
well-known radio-IR correla tion seen for brighter galaxies 
l|Bellll2003l : iBovle et allhoojj ) . 



Space Observatory {ISO) allowed fainter and higher red- 
shift galaxies to be observed at mid-infrared (MIR) wave- 
lengths. These ISO observations showed that the MIR emis- 
sion from galaxies is loosely correlated with radio emission 
acros s a wide range of redshifts, tentatively extending to 
z~4 l|Cohen et aZj|2000l : lElbaz ef fflZj|2002llGarrett 112003 ). 

More recently the launch of the Spitzer Space Tele- 
scope in August 2003 has greatly increased the sensitivity of 
MIR observations and hence our ability to study the MIR- 
radio correlation. Early results, such as from the Spitzer 
First Look Survey (FLS), have confirmed that the MIR- 
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iBel ll2003h argue that while the IR emission from lumi- 
nous galaxies will trace the majority of the star-formation 
in these sources, in low luminosity galaxies the IR emis- 
sion will be less luminous than expected considering the 
rate of star-formation within the source (i.e. the IR emis- 
sion will not fully trace the star-formation). In this scenario 
the reduced efficiency of IR production relative to the source 
star-formation rate (SFR) would be the result of inherently 
lower dust opacities in lower luminosity sources and conse- 
quently less efficient reprocessing of UV photons from hot 
young stars into IR emission. The simple consequence of 
this is that at lower luminosities the near linear radio-IR 
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Figure 1. Left-hand panel: Radio 1.4 GHz versus the MIR 24 fim flux density of all 377 individual sources (small triangle), and median 
radio flux density logarithmically binned by their 24 /an flux density (filled circles) within the 8'5x8'5 field. Right-hand panel: Control 
plot of 1.4 GHz radio flux densities plotted against the 24 /im source flux densities of the sample. The radio flux densities of the control 
sample have derived in exactly the same manner as the flux densities plotted in the left-hand, panel but have been measured at randomly 
assigned sky positions in the 8'.5x8'.5 radio image where no known sources at any wavelength are located. 



correlation L ra dio c<Lj R , with 7 > 1 (e.g ICox et al .1 1 19881 : 
IPrice fc Duric Ill992t) will be deviated from. Of course such 
an assertion is dependent upon the radio emission providing 
a reliable tracer of star-formation at low luminosities which 
may be equally invalid. 



2 DATA AND ANALYSIS 



2.1 MERLIN+VLA observations 



Recently iBovle et al\ (|2007n have presented a statisti- 
cal analysis of Australia Telescope Compact Array (ATCA) 
20 cm observations of the 24 fim sources within two regions 
(the Chandra Deep Field South (CDFS) and the European 
Large Area ISO Survey SI (ELAI S)) of the Spitzer Wide 
Field Survey (SWIRE). In this work lBovle et al\ (|2007f ) have 
co-added sensitive (rms~30 /^Jy) radio data at the locations 
of several thousand 24 /j,m sources. Using this method they 
have statistically detected the microJy radio counterparts of 
faint 24 /j,m sources. At low flux densities (S24nm = 100 /iJy) 
they confirm the IR-radio correlation but find it to have 
a lower coefficient (Si.4GHz = 0.039 S24 ^m) than had previ- 
ously been reported at higher flux densities. This coefficient 
is significantly different from results p reviously derived from 
detections of individ ual objects (e.g. Appleton et aZ.I [20041 ') 
and is speculated bv lBovle et al\ ||2007h to be the result of 
a change in the slope of the radio-IR correlation at low flux 
densities. 

In this paper, we utilise very deep, high resolution 20 cm 
observations of the Hubble Deep Field North and surround- 
ing a rea made using MERLIN and the VLA |Muxlow et al\ 
120051 ) in combination with publicly available 24 /jm Spitzer 
source catalogues from GOODS to study the MIR-Radio 
correlation for microjansky radio sources. This study ex- 
tends the flux density limits of the radio-IR correlation by 
more than an order of magnitude for individual sources and 
overlaps the flux de nsity regime studie d using statistical 
stacking methods by IBovle et al\ (|2007l ). Additionally we 
em ploy statisti c al sta cking methods, similar to those used 
bv lBovle et~ai\ (|2007r i, to extend the correlation further to 
still lower flux densities. 



We adopt H =75kms _1 Mpc~\ fi n 
throughout this paper. 



: 0.3 and f2 A =0.7 



Extremely deep radio observations of the HDF-N region 
were made in 1996-97 at 1.4 GHz using both MERLIN and 
the VLA. Th e se ob s ervations were i nitiall y pre sented in 
iMuxlow etd\ (120051 ). iRichards et~ai\ |l998T ) and iRichardsl 
l|200d ). Theresults from the combined 18 day MERLIN 
and 42 hr VLA ob servations are described in detail in 
IMuxlow et'oil (|2005l ). The combined MERLIN+VLA image 
has an rms noise level of 3.3 /iJy per 0'.'2 circular beam mak- 
ing it amongst the most sensitive, high-resolution radio im- 
ages made to date. 

Using the same methods as described in lMuxlow et al] 

({2JD05J) these combined MERLIN and VLA observations 
have been used to image the entire unaberrated field of 
view, 8.5x8.5 arcmin 2 in size, centred upon the original 
MERLIN pointing position (a = 12 h 36 m 49=4000, 8 = 
+62° 12' 58'.' 000 (J2000)f]. This image has an rms noise level 
of 3.6 ^iJy beam -1 and has been convolved with a 0'.'4 cir- 
cular beam. These observations have been shown to align 
with the International Coordinate Referenc e Frame (ICRF) 
to better than 15 mas (|Muxlow et alj|2005h . 



2.2 GOODS-N Spitzer 24 (im observations 

As part of the GOODS enhanced data releas^E] (DR1+ 
February 2005) a catalogue of Spitzer 24 /jm source posi- 
tions and flux densities for the GOODS-N field were re- 
leased. This source catalogue is limited to flux densities 
>80/iJy providing a highly complete and reliable sample. 
At the time of writing this 24 /im source catalogue repre- 



1 Public access to these radio data will be made available via 
an on-demand radio imaging tool developed by the authors using 
Virtual Observ atory tools. A detailed explanation of this service 
can be found in Richards et al. (2007). 

2 http:/ /www. stsci.edu/science/goods/DataProducts/ 



An evolution of the IR- Radio correlation? 3 




100 1000 10000 

S 2 4„m (MO 



Figure 2. Radio 1.4 GHz versus the 24 /an flux density. Sources 
from the 8'.5x8'.5 HDF-N field are plotted individually (small 
black triangles). The median radio flux density logarithmically 
binned by 24 fim flux densities are plotted as large filled red cir- 
cles. The solid red line represents the best-fit line to the binned 
HDF-N data alone. Sources within the CDFS-SWIRE field de- 
tected at both 24 (im and 1.4 GHz from Norris et al. (2006) 
are plotted as either green open stars (identified as AGN) or 
green triangles (not identified as AGN). All sources detected at 
both 1.4 GHz and 24 /im in the Spitzer First Look Survey (FLS) 
with source positio n separations of <1"5 are plotted in orange 
llFadda et adl2006l) . Note the quantisation of the SWIRE and 
FLS points, in this and subsequent plots is a result of the ac- 
curacy of the flux densities tabulated in the literature. The blue 
pluses and fitted lines in the low portion of the plot show the 
IR-radio correlation derived from stacking radio emission at the 
positions of 24 (im sources in the CDFS and ELAIS field by Boyle 
et al. (2007). A brief summary of the basic characteristics of each 
of the data sets plotted included in TablcfT] 



sents the most complete and accurate mid-infrared source 
list publicly available for the GOODS-N/HDF-N field. 

All 24 /im sources which are detected optically in 
GOODS HST ACS images show an accurate astrometric 
alignment with their optical counterparts implying that the 
astrometry between these two data-sets and their subse- 
quent catalogues is self-consistent. However, a comparison 
of the astrometric alignment of the positions of sources cat- 
alogued by GOODS derived from their HST ACS images 
(Richards et al. 2007; Muxlow et al. in prep) shows there 
to be a systematic offset in declination of — 0'.'342 from the 
radio reference frame. This linear declination correction, al- 
though small relative to the Spitzer resolution at 24 /im is 
significant when compared to these high resolution radio 
data. This linear correction has been applied to the Spitzer 
source positions prior to all comparisons between the two 
data sets. 



2.3 Measurement of source radio flux densities 

The field covered by these GOODS-N Spitzer 24 pun observa- 
tions contains 1199 24- fim sources identified with flux den- 
sities >80piJy. The 8.5x8.5 arcmin 2 radio field contains 377 
of these Spitzer sources with 24 fim flux densities ranging 
from 80.1 to 1480 fiJy. 

Radio emission with the deep MERLIN+VLA data has 
been measured at the corrected position of each individual 
24 pun source within the radio image. The radio flux den- 
sity was measured within a series of concentric rings with a 
maximum radius of 4arcsec centred on each 24 pm source. 
Statistical analysis of these measurements and a sample of 
stronger radio sources within the same field shows that the 
radio flux density increases within progressively larger an- 
nuli out to a radius of l'.'5. For stronger (> 5a) radio sources 
the total flux density recovered using this method within a 
radius of l'/5 is between 90% and 100% of the source's total 
flux density as measured by other means, such as Gaussian 
fitting. Additionally the average (either mean or median) ra- 
dial profile of all of the 24 /im sources shows that almost all 
of t he radio flux density is recovered within radii less than 
l'/5 (IBeswick et al. I l2006l ; Beswick et al. in prep) . Through- 
out the rest of this paper the radio flux density recovered 
within a radius of l'.'5 of each IR Spitzer source will be as- 
sumed to be equal to the total radio source flux density. 



As a confirmation of this method several thousand ran- 
dom sky positions within the radio field were selected and 
cross-referenced to exclude any which were at the positions 
of known sources at any wavelengths. At each of these blank 
positions the radio flux density was measured using an iden- 
tical method as outlined above. The flux densities recorded 
within each of these annuli showed no positive bias within 
any ring and followed the same Gaussian distribution as de- 
rived from the pixel-by-pixel noise statistics of the whole 
image once all radio sources greater than 40 /iJy have been 
subtracted. As further test of these methods 90 fake, faint 
(< 3<r) and partially-resolved radio sources were injected 
into the u-v data. These data were then Fourier inverted and 
cleaned in an identical manner to the real data and then the 
flux densities of these fake sources were extracted using the 
methods outlined above. Whilst, as expected, these sources 
were not detected individually their radio emission, when 
averaged together, was consistent with the average flux den- 
sity of the population of fake sources injected into the data. 

A more detailed description of this method and the ra- 
dio characteristics of this sample of 24 pun Spitzer source, 
including source catalogues and size information, will be pre- 
sented in a forthcoming paper. 



3 RESULTS 

Of the 377 24 pirn Spitzer sources within the radio field 303 
were found, by this method, to have total radio flux densities 
greater than 3 times the local map rms. Many of these source 
have peak flux densities below 3-er and hence were not se- 
lected using the radi o-only detection methods employed by 
iMuxlow et ail 12005). The 1.4 GHz versus 24 pun flux densi- 
ties of all 377 of these source are plotted in left-hand panel 
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Table 1. Summary of deep 1.4 GHz radio and 24 fira surveys 


refereed to in this paper. 




Survey 


Radio 


Instrument 


Radio senitivity 


Angular 


Reference 




Survey Area 




(/ijy beam ) 


Resolution 




HDF-N / GOODS-N 


72.25 arcmin 2 


VLA+MERLIN 


3.3 


(?.'2x(f.'2 


Muxlow et al. 2005, this paper 


ATLAS (CDFS/ELAIS) 


3.7 deg 2 


ATCA 


20^60 


11" x5" 


Norris et al. 2006 


Spitzer FLS 


5 deg 2 


VLA 


23 


5" 


Condon et al. 2003, Fadda et al. 2006 
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Figure 3. Mean (left) and median (right) images of the 1.4 GHz radio emission for all sources within the six faintest 24 /im flux density 
logarithmic bins plotted in Fig. 1 in descending flux density order from top-left to bottom-right. The range of 24 fim flux density over 
which each image has been stacked is shown at the top of individual panels. Each image is contoured with levels of —2, —1.414, —1,1, 
1.414, 2, 2.828, 4, 5.657, 8, 11.31, 16, 22.63 and 32 times 3x(3.3/v / N) /iJybm^ 1 , where N equals the number of 24 jim source positions 
averaged in the map. The peak flux density, lowest plotted contour and number of IR sources which have been averaged over (N) is listed 
at the bottom of each image panel. 



of Fig[T] The control sample, plotted in the right-hand panel 
of Fig[T] is derived using an identical radio flux extraction 
methods and from the same radio data but at randomly 
assigned locations coincident with no known source posi- 
tion. The flux densities of these sources are also plotted in 
Fig(2] along with sources d etected in the shallo wer surveys of 
the CDFS-SWIRE field bv lNorris et aZI d2006f) using ATCA , 
and the Spitzer FLS using the VLA bv lFadda et all (|200fJ) . 
Sources from the so uthern CDF S -SWI RE field which have 
been categorised bv lNorris etUll (|2006j ) as containing AGN 
(plotted as stars) are generally situated toward the upper 



half of this figure and show an expected radio excess when 
compared to their 24 /xm flux densitie s . Sou rce which were 
not identified as AGN bv lNorris et all (|2006l ) are plotted as 
triangles in this figure. It should be noted that these sources 
have not been classified as containing an AGN but may in 
fact be either AGN, star-forming galaxies or a combination 
of the two. 

Overlaid in blue in the lower half of Fig. [2] are the flux 
densities of groups of median stac ked infrared source s from 
the ELAIS and CDFS field from iBovle et all <|2007l ) along 
with fits to these data. Also overlaid, as filled red circles, in 
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this figure are the median values of all of the HDF-N sources 
logarithmically binned by their 24 /j,m flux densities. There is 
a clear discrepancy between the averaged radio flux densities 
derived from these HDF-N observations and the flux densi - 
ties from the ELAIS and CDFS fields l|Bovle et all 120071 ') 
which will be discussed further in the following sections. 

Using an analysis method similar to lBovle et all (|2007h 
the mean and median images of the radio emission from 
HDF-N field of the 6 lowest 24 fim flux density logarith- 
mic bins (as plotted in Figs.[T] and [2| have been derived 
and are shown in Fig. [3] Each of these images is the sta- 
tistically combined radio emission from the location of sev- 
eral tens of Spitzer 24 /im sources and has been contoured 
at multiples of three times 3.3 /iJy beam -1 divided by %/N 
where N equals the number of Spitzer source positions that 
have been stacked together. As can be seen in these images 
the off-source noise levels achieved approaches the value ex- 
pected when co-adding multiple images with near-Gaussian 
noise properties. The co-added image rms achieved in the 
faintest 24 /im flux density bin (80.0 to 100.25 /iJy) is 0.45 
and 0.56 /iJy beam -1 in the mean and median co-added im- 
ages respectively. The flux densities measured directly from 
these stacked radio images are consistent with the statistical 
mean and median flux densities within the equivalent bins 
derived from the individual source flux extraction method 
described in Section 2.3 (Table[2]) and plotted in Figs. [1] and 
[3 In each stacked image the largest angular size of the re- 
gions showing elevated radio emission within Fig. [3] is en- 
closed within the l'.'5 radius circular area from which the 
individual source flux densities have been extracted. 

The sizes of the composite radio sources in these co- 
added images (Fig.[3]) represent a combination of the true 
radio source size, any offsets between the IR and radio emit- 
ting regions within the sources and any random errors in 
the Spitzer source positions, all of which will be convolved 
with the synthesised beam applied to these radio data (0'.'4 
circular beam in this case). No systematic offset is seen be- 
tween the emission in these composite radio images from 
their nominal image centres implying that no significant 
non-random errors between the astrometric alignment of 
these data-sets remain. Further analysis of the radio source 
sizes and structures will be the subject of a future paper. 

In Fig. [3] the mean (left) images of the radio emis- 
sion for the two highest 24 /im flux density bins are dom- 
inated by emission from a few bright sources. Several of 
these brighter radio sources can be seen in Fig. [T] scattered 
above Si.4GHz =200 fiJy. The presence of these brighter ra- 
dio sources and the low numbers of images stacked together 
is clearly evident when the mean and median averaged im- 
ages are compared. In particular in the second brightest bin 
(197.27 to 247.27 ^Jy) two mjy radio sources are present and 
offset from their equivalent 24 /im position. It is likely that 
these bright radio sources, which deviate from the radio-IR 
correlation, may contain radio AGN which would feasibly be 
spatially offset from any ongoing star-formation and 24 /im 
position. 

An alternative method of projecting the correlation 
between these 24 and 1.4 GHz data is as a function 
of the commonly defined parameter q24, where q-24 = 
log(S24/j m /Si.4GHz). The mean and median value of q24 for 
all of the HDF-N sources presented in this study is 0.52 and 
0.48 respectively, compared to the value of 0.69±0.36 and 



0.73 derived from non-AGN sources catalogued in observa- 
tions of the SWIRE field by iNorris et all j2006l) . A similar 
study of brighter sources in the Spitzer FLS (jAppleton et all 
I2004T ) show the mean value to be (q 24 =0.84±0.28) slightly 
higher than average value derived here. The distribution of 
q24 values for all sources in this HDF-N study is shown in 
Fig.H 

In Figs.[5][6][7]and[8]the parameter q 2 4 is plotted against 
S24/jm, Si.4GHz and redshift. The monochromatic value of 
q24 represents the slope of the IR-radio correlation and the 
size of dispersion in q24 is related to the strength of the 
correlation. In each of these figures no fc-correction has been 
applied to either the radio or IR luminosities since we have 
incomplete redshifts for the sample, and inadequate spectral 
information is known about the majority of the individual 
sources in either the IR or radio bands. A similar study by 
lAppleton et all (|2004T ) has shown that fc-corrections made to 
individual sources result in only a minimal systematic effect 
on the average value of q24 below redshift ~1. 

In Fig.[5]the values of q24 for sources fr om these HDF-N 
data, the SWIRE field (jNorris et aZj|2006h and the Spitzer 
FLS field (|Fadda et a/JbOOrJ ) are plotted against 24 /im flux 
density. In the uppe r panel both the A GN identified and 
non-AGN sources of INorris et all (|200ot ) are shown whilst 
in the lower panel the AGN sources are excluded for clar- 
ity. In the upper panel the AGN sources from the SWIRE 
field, as expected, deviate from the correlation generally 
showing an excess of 1.4 GHz radio emission compared to 
their 24 fim emission. The HDF-N field covers a consider- 
ably smaller area and was originally chosen to have few 
bright radio- AGN, which is confirmed by the minimal num- 
ber of HDF-N sources displaying an excess of 1.4 GHz radio 
emission. Amongst the HDF- N sources and the non-AGN 
sources from the SWIRE field (jNorris et aZj|2006l . and Fig.[5] 
here) q24 shows a small trend toward lower values at smaller 
24 fim flux densities. The SWIRE sources plotted here have 
been selected to have both radio and IR flux density above 
a certain threshold. The lowest radio flux density of sour ces 
within the SWIRE field catalogue (|Norris et all l200fj ) is 
100 piJy. A line of constant radio flux density (100 /iJy, the 
lowest individual source detections in the SWIRE catelogue) 
is overlaid on Fig.0 above which sources from the SWIRE 
study are excluded. The HDF-N radio data plotted here has 
a sensitivity limit ~10 times lower than that of the SWIRE 
data, thus potentially only excluding sources in the extreme 
top left-hand part of this plot. Furthermore, the sample of 
sources plotted here has been subject to no radio detec- 
tion threshold limit hence removing this potential bias com- 
pletely. 

The q 2 4 values for the HDF-N data and the sources from 
INorris eta l. (2006) binned as a function of S24 M m, along 
with the ratio of 824^111 and Si.4GHz are shown in Fig[6] 
Especially within the binned points, these two plots show a 
tentative trend for low q24 and C S24 with declining 24 urn 

»1.4GH. r 

flux density. The gray areas and dotted line in these two pan- 
els show the value of t he q24 as derived from the Spitzer FLS 
l|Appleton et all 12003 ). The HDF -N observations presen ted 
here are consistent with those of lAppleton et all |2004l ') at 
high 24pm flux densities. 

The effect of the 24 /im sensitivity cut-off (80.1 A*Jy) in 
the GOODS - N dat a and the equivalent limit in the data of 
INorris et a l. (2006) will, however, positively bias measure- 
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Table 2. Flux measurements from stacked images and binned data. The flux density in the stacked images has been measured by fitting 
a single Gaussian to the radio emission. The errors attached to measurements made from the stacked images are the formal errors of 
these Gaussian fits, rather than the image rms which in each case is comparable to a. Sources within 3arcsec of the edge of the radio 
field have been excluded from the stacked images, hence slightly reducing the number of sources combined in each image relative to those 
included in the binned data. 



Binned data Stacked images 







Median 


Mean 






Median 


Mean 


S24 |im 


Number 


Sl.4GHz 


Sl.4GHz 


CT 1.4 GHz 


Number 


Si. 4 GHz 


Sl.4GHz 


(MJy) 




(pJy) 


(A*Jy) 


(A«Jy) 




G"Jy) 


(MJy) 


88.8 


65 


28.9 


22.6 


0.68 


65 


26.7±8 


27.8±6 


112.5 


72 


34.1 


32.0 


0.62 


70 


34.0±7 
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ments of q24 at low radio flux densities. This is clearly shown 
in Fig. [7] The lower flux density limit imposed by the Spitzer 
sensitivity results in an exclusion of sources in the lower por- 
tion of Fig. [7] This consequently biases the determination of 
the value of q24 for any complete sample of sources with low 
radio flux densities. 

Of the 377 Spitzer source within the 8'.5x8'.5 HDF-N ra- 
dio field, 259 of these galaxies have published spectroscopic 
or photometric redshifts. The q24 values for this subset of 
sources, along with 50 non-AGN SWIRE sources, are plot- 
ted with respect to their redshift in Fig[S] Within these two 
data-sets q24 is seen to slightly reduce at higher redshifts 
although this effect is small and significantly less than the 
scatter. These values of q24 have not been ^-corrected. The 
application of a ^-correction to these data will result in a 
small increase in q24 which will increase as a function of z. It 
is interesting to also note that recent studies o f discrete areas 
within 4 very nearby star-forming galaxies l|Murphv et al\ 
2006a) have show n q24 to vary across the extent of individual 
sources. However, iMurphv et all(|2006al ) derive a mean value 
of q24~l when integrating over entire galaxies in their sam- 
ple, at z«0, which is consistent with many of the source in 
both the SWIRE and HDF-N samples but somewhat higher 
than the average value derived at higher redshifts. 

Using the available redshifts the observed luminosity 
has been calculated and is plotted in Fig. including both 
the sources in the HDF-N region and the non-AGN sources 
in the SWIRE field. As can be seen in Fig. [9] the 24 /xm-to- 
1.4 GHz correlation, for individual sources, extends down to 
Li.4ghz » 10 21 WHz- 1 (L24 fim « lO^WHz- 1 ). 



4 DISCUSSION 

4.1 Radio emission from faint 24 fim sources:- 
Extending the IR-Radio correlation 

The primary results of this work demonstrate that this sam- 
ple of individual 24 /im selected sources follow the MIR-radio 
correlation down to radio flux densities of a few microjy. 
These results are consistent with shallower independent ra- 
dio surveys (such as the ATCA observations of the SWIRE 
field) and seamlessly extend the IR-radio correlation down 
to lower flux densities (see Fig. [2}. 




Figure 4. Histogram of the distribution of q24 for the 377 sources 
within the 8'. 5x8'. 5 field centred upon the HDF-N. 



The majority of the IR selected sources within the HDF- 
N region have radio flux densities less than 100 fiJy, with 
only a few sources which show a significant excess of 1.4 GHz 
flux density, when compared to the IR-radio correlation. 
This implies that the vast majority of these IR-selected mi- 
croJy radio sources are primarily driven by star-formation 
with little significant conta mination from AGN. T his is con- 
sistent with the results of iMuxlow et al\ (2005) who find 
that below 100 /iJy greater than 70 percent of radio sources 
are starburst systems. However, the data presented here are 
based on an IR selected sample rather than a radio selected 
sample. This will inherently result in the sample being domi- 
nated by star-forming systems but will not exclude the pres- 
ence of a population of weak microjy AGN sources. 

At the positions of the majority of the lowest 24 pirn flux 
densities sources the peak radio emission is not significantly 
greater than the image rms. As a consequence most of these 
source s were not formally identified in the previous radio 
study (IMuxlow et qZ.ll2005h . In many case, especially at low 
24 fim flux densities, the extracted total radio flux density 
of individual sources is close to the radio sensitivity of these 
data (Fig[T]). When directly compared with the control sam- 
ple (Figfl] Right), which has been compiled by extracting the 
radio flux density from blank areas within the same data us- 
ing identical methods, a significant correlation between the 
strength of the total radio and 24 /im emission is evident. 
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Figure 5. The 24 fim flux density versus q24. All 377 source from 
the 8'.5x8'.5 HDF-N field are plotted as individual sources (small 
black triangle). In the upper panel sources detected at both 24 (im 
and 1.4 GHz in the SWIRE-CDFS field from Norris et al. (2006) 
are plotted in green (small triangles are sources not identified as 
AGN and open stars have been identified as AGN). Sources from 
the Spitzer FLS are plotted in orange from Fadda et al. (2006). In 
the lower panel all of the Spitzer identified sources in the HDF- 
N and the only non-AGN sources from Norris et al. (2006) are 
plotted. Additionally plotted as 'pluses' and 'crosses' in the lower 
panel are the values of q24 as derived from the stacking analysis 
by Boyle et al. (2007) of the SWIRE-CDFS and ELAIS fields 
respectively. The solid green line overlaid on both panels is a line 
of constant radio flux density of 100 /iJy, the lowest radio flux 
density of sources in the SWIRE sample plotted here. 



However, the radio sensitivity limit of these data precludes 
the detailed discussion of individual low flux density sources. 

4-1.1 The average radio emission from the faintest IR 
sources 

The radio emission from many of the individual low flux 
density 24 /im sources is too faint to be imaged with high 
signal-to-noise in these radio data. However, it is possible to 
characterise their average radio emission by either binning 
the measured radio excess at their locations or by forming 
composite radio images from many sources. With the rela- 
tively limited sample of IR sources which are co-spatial with 
these radio observations (377 in total) it is only statistically 
viable to use either binning or image stacking methods at 
the lowest IR flux densities where the number of sources 
that can be combined becomes large. 

Overlaid in Figs. [T] and [2] are the median flux densities 
of IR sources in the HDF-N region imaged here. These bins 
have been logarithmically sampled as a function of 24 /im 
flux density and include all sources within these bins re- 
gardless of their measured radio flux density. Whilst at the 



Figure 6. In the upper panel the flux density ratio 
(S24 /jm/Si.4 GHz) versus 24 /im flux density is shown. The in- 
dividual sources from the HDF-N field are plotted as small black 
triangles, the median values of these HDF-N source binned as a 
function of S24 as filled circles, green stars show the median 
binned values for sources listed as non-AGN within the CDFS- 
SWIRE sample of Norris et al. (2006) and the blue 'crosses' and 
'pluses' show the flux density ratios derived from the stacking 
analysis of the SWIRE-CDFS and ELAIS fields respectively from 
Boyle et al. (2007). The overlaid black dotted line is the mean 
value for (S24 ^m/Si.4 GHz) derived by Appleton et al. (2004). 
This line is equivalent to q24 =0.84±0.28 with the gray filled box 
representing the area enclosed by these errors and the flux den- 
sity range investigated by Appleton et al. (2004). The values of 
q24 against 24 fira flux density are plotted in the lower panel. 
The symbols within this plot are identical to those in the upper 
panel, individual HDF-N sources are not included for clarity. The 
additional diagonal solid green line represents a line of constant 
radio flux density of 100 /tJy, the lowest flux density of sources 
in the SWIRE sample plotted here (sources above this line are 
excluded by this limit from the SWIRE sample). This flux den- 
sity limit will significantly effect the values of the binned SWIRE 
data points (green crosses) negatively biasing the four lowest flux 
density bins. This bias only effects the SWIRE sample. 



lowest emission levels the radio flux densities of individ- 
ual sources will be subject to moderate or large errors the 
binned ensemble of these points will robustly represent the 
flux distribution of sample as a whole. Using the radio flux 
densities extracted at the positions of the individual 24 /im 
sources, significant radio emission is detected in every aver- 
aged 24 fim flux density bin. 

To further test the viability of measuring the statisti- 
cal radio flux density from sources in this sample below the 
formal radio detection limit of these MERLIN+VLA data, 
the images have been stacked at the positions of multiple 
IR sources, binned as a function of their 24 /im flux density 
(Fig.[3)l . These postage-stamp images have been averaged on 
a pixel-by-pixel basis, where each pixel in the stacked image 
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Figure 7. The 1.4 GHz flux density versus q24. The samples 
plotted and symbols used are the same as those plotted in 
Fig[5] In both plots the solid line represents a line of constant 
S24/im=80.1 fijy which is equal to the detection threshold of the 
GOODS-N Spitzer sample used. 



is either the mean or median value of the pixels within the 
individual 24 /im source image stacks. Each of the 6 low- 
est bins, for both the mean and median averaged images, 
show significant levels of radio emission. The 1.4 GHz flux 
density in these images, extracted via Gaussian fitting, is 
consistent with the binned level of emission derived from 
measuring individual sources, demonstrating the validity of 
either method (Tabled . 

The sub-arcsecond angular resolution of these MER- 
LIN-!- VLA observations, in conjunction with their relatively 
good u-v coverage, allow the average radio structures of faint 
IR sources to be investigated. The median and mean images 
in Fig.|3j especially in the two highest flux density bins show 
some distinct differences. These differences arise due to in- 
clusion of a few bright radio sources which dominate the 
average emission included in these bins. In particular in the 
second highest flux density mean image several 'bright' ra- 
dio sources are offset from the map centre. The positional 
offset of radio emission from these sources is real and in- 
dicates that, at least in these cases, that the origins of the 
radio and IR emission are different. 

The Gaussian fitted sizes of the radio emission in the 
stacked images (Fig[3} provide an upper limit on the aver- 
age size of the radio counterparts of these faint IR sources. 
The largest angular sizes of the radio emission in the median 
stacked images created from this sample range between l'/4 
and 2". This is approximately equivalent to a linear size of 
lOkpc at redshifts beyond 1. These upper limits on the ra- 
dio source sizes are consistent with radio emission on galac- 
tic and sub-galactic scales and originating within kpc-scale 
starburst systems. A more detailed analysis of the sizes and 



Figure 8. The q24 versus redshift. All 259 source from the 
8'.5x8'.5 HDF-N field are plotted as individual sources (black 
triangle) with known redshifts. Fifty sources from the CDFS- 
SWIRE field detected at both 24 /im and 1.4 GHz, and redshift 
information from Norris et al. (2006) are plotted as green tri- 
angle. Only sources not identified as AGN by Norris et al. 2006 
are plotted. No k-correction has been applied to the data points 
plotted. 



structures of the radio emission from both these composite 
images (Fig.[3jl and the individual sources themselves is be- 
yond the scope of this paper and will be presented in a latter 
publication. 



4-1-2 Comparison with other results 

Two directly comparable studies of the 1.4 GHz radio and 
24 /im MIR emission from sources in t he Spitzer FLS 



Applcton et al. 



Appleton et al. 



and t he S WIRE fields have b een made by 
(|2004j) and iBovle et all (j2007l ) respectively. 
(2004 ) compa r ed ob servations of the Spitzer FLS made by 
ICondon et all (|2003l ) using the VLA in its B-configuration 
(5"restoring beam) with 24 /im and 70 /im Spitzer observa- 
tions. These observations had a detection threshold of 90 ^tJy 
and 500/xJy at 1 .4 GHz and 24 pm respectively. Whereas 
IBovle et all (|2007l ) co-added deep (a ~30/jJy/bm) 1.4 GHz 
ATCA observations of the southern SWIRE field at the posi- 
tions of many thousands of Spitzer 24 fim sources, regardless 
of their detection at radio wavelengths, in order to statisti- 
cally detect radio emission at the few fijy level. As a compar- 
ison to the HDF-N obs ervations presented here, data points 
from the SWIRE f ield (|Norris et ffll.ll2006f ). the Spitzer FLS 
dFadda et aZj|2006h and the co-added flux densities derived 
bv lBovle et all (|2007i ) are pl otted in FigsJl G p Ela nd [71 

Using the Spitzer FLS I Appleton et all l|2004l ) studied 
MIR to radio correlation for 508 sources with redshift in- 
formation, confirmin g results from simila r previous surveys. 
Following this work lFadda et al] (|2006t ) have published a 
more comprehensive catalogue of ~17,000 24 fim sources 
from the Spitzer FLS of which 2415 have radio counter- 
parts within 4". Just over 2000 sources with radio counter- 
parts within l'.'5 are plotted in Figs. [2] [5] and [7] In general 
the sources plotted follow the MIR to radio correlation but 
these data contain many outliers with elevated radio emis- 
sion, characteristic of galaxies containing AGN. 

The study of lAppleton et all (|2004l ) included sources 
with S24/J™ ~0.5mJy and Si.4GHz >90 fiJy. The majority of 
this sample have redshi fts of ~1. The mean (no n-fc-corrected) 
value of q 2 4 derived by lAppleton et all l|2004h is 0.84±0.28, 
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which is shown in Fig[6] The non-AGN sources detected 
at both 1 4 GHz and 24 um w ithin the SWIRE field cata- 
logued bv iNorris et all (|2006h are also plotted in Figs. [2] [5] 
and [7] These sources are limited to ^0.1 mjy and >170/iJy 
at 1.4 GHz and 24 pm respectively. Using only the sources 
within the SWIRE field identified as not being AGN by 
INorris et all l|2006l ) the mean (non-fc-corrected) values of q24 
is 0.69±0.39. The radio flux density threshold of both of 
these surveys overlaps with the brighter radio sources in 
this MERLIN+VLA study. 

These deep HDF-N MERLIN+VLA observations trace 
individual sources with considerably fainter flux density 
than those observed in the SWIRE and Spitzer FLS sur- 
veys. These observations (Figs. [2] [5] and [7| are consistent 
with the extrapolation to lower flux densit i es of results 
found previously (e.g. Appleton et all |2004| ; INorris et all 



l200fj ; iFadda et aZ"Tl2006t ). However both the SWIRE and 
Spitzer FLS surveys which sample sources with somewhat 
higher flux densities sources have mean values of q24 which 
are consistent but slightly higher than the value of 0.52±0.37 
derived from these HDF-N data. 

Averaging together the radio emission from multi- 
ple 2 4 lira Spitzer sourc e within the CDFS and ELAIS 
fields iBovle et all (|2007m have traced the radio emission 
against MIR sources with flux densities comparable to 
those detected in t hese more sensitive HDF-N observations. 
IBovle et all l|2007l ) have exploited the large number of MIR 
sources detected in these fields to average together multiple 
radio images made at the MIR source positions and statisti- 
cally detect radio emission significantly below their radio de- 
tection threshold. This method is directly comparable to the 
methods applied to these radio observations of the HDF-N 
and is identical to the methods used to produce the stacked 
images presented in Fig. [3] and Table[2] The me dian stacked 
radio flux densities of MIR sources derived by IBovle et all 
(2007) show a significantly reduced Si.4GHz relative to S24fim 
(i.e. higher q24 value) compared with the distri bution of in- 
dividual sources detected in the Spitzer FLS (I Fadda et all 
l200fj: I Appleton et qZ.ll2004h . the catalogued SWIRE sources 
|Norris et aZ.I l2006h or the HDF-N results presented here. 
Whilst ther e are s everal differences between the studies of 
IBovle et all (|2007l ) and those presented here, such as the 
sensitivity (rms ~30 /iJy versus ~3 /^Jy) and angular resolu- 
tion (6"versus 0'.'4) of the two surveys, it remains hard to 
simp ly reconcile these results. 

IBovle et all (|2007h recognised the inconsistency of their 
calculated values of q24(=1.39) when compared with 0.84 
derived bv I Appleton et all (|2004T l and have extensively sim- 
ulated and tested the correctness of thei r analysis meth- 
ods. Since the results of lBovle et all l|2007l ) are based on an 
infrared-selected sample with no radio detection limit, this 
discrepancy could possibly be explained by the presence of 
a previously unknown population of faint radio sources with 
a low radio to infrared flux densi ty ratio. However the val- 
ues derived by IBovle et all l|2007l ) are not confirmed in the 
study presented here. 

A possible alternative explanation may be that the 
ATCA observations are underestimating the radio flux den- 
sity of sources at the lowest flux densities when compared 
to these other surveys which have used radio data ob- 
tained using the VLA and/or MERLIN. It is interesting 
to note that an earlier deep survey using the ATCA at 



1400 MHz (ATESP, A ustralia Telescope ESO Slice Project) 
IPrandoni etall (|2000h found, upon comparison with NVSS 
flux densities for the same field, that whilst for bright sources 
the ATCA and NVSS flux density values were consistent, 
for faint sources the ATCA results underestimated the ra- 
dio flux densities by up to a factor of 2 compared with the 
VLA N VSS. This underestima tion of radio flux density was 
scon bv lPrandoni et all (2000) to increase with diminishing 
radio brightness. Whilst it remains unclear as to the cause, 
or reality, of this effect, if applicable to the ELAIS and CDFS 
fields, it could res ult in the elevated values of q24 seen by 
IBovle et all l|2007l ) at very low radio flux densities. 

4.2 An evolution in the IR-radio correlation? 



It has been sug gested, both theo retically (Bell 2003j) and 
observationally (|Bovle et aLll2007f ). that the IR-radio corre- 
lation may evolve at lower flux densities. One of the aims of 
this study is to investigate this correlation for very low flux 
density radio sources. 

Derived from the HDF-N data alone the average value of 
q24 shows a small trend to reduce as a function of lower val- 
ues of S24,jm (Fig.[5]&[6]). This can be further demonstrated 
by splitting these data into two approximately equal sam- 
ples, above and below S24,jm=150/iJy. The mean and me- 
dian values of q24 for the 187 sources below 150 pJy are 0.45 
and 0.36 respectively, compared with 0.56 and 0.54 derived 
from the 190 sources with S24 M m >150/xJy. These HDF-N 
results along with res u lts from samples o f brighter galaxies 
l|Appleton et all |2004| ; INorris et all |2006|) are broadly con- 
sistent with average q24 reducing as a function of decreasing 
S 2 4 Mm (see Fig.[5]&[l. 

There are several contributing factors which will affect 
these results and those from previous studies. These factors 
include biases introduced by flux density thresholds of obser- 
vations, the application or non-application of fc-corrections 
and possible AGN contaminat ion of sample ga l axies. 

Previous studies, such as lAppleton et all f2004), have 
used samples of sources detected at both radio and infrared 
wavelengths. The implications of a radio detection limit on 
these results will be to exclude sources with high q24 at low 
S24/jm. The line overlaid on Fig. [S] represents the upper limit 
imposed on the ATCA observations of the S WIRE survey 



by th e radio detection threshold of these data (INorris et al. 



2006). In a manner similar to that used by IBovle et al. 



( 2007) no radio detection limit has been applied during anal- 
ysis of these MERLIN+VLA data with the radio flux density 
being recorded for all MIR sources regardless of significance 
of their radio detection. This method will result in errors in 
the measured radio flux densities for very weak individual 
radio sources but will not directly bias against larger values 
of q24. However when many sources are statistically aver- 
aged together the mean and median value of detected radio 
emission and hence the derived value of q24 will be robust 
even for source considerably fainter than the radio detection 
threshold of these data. 

No direct fc-correction has been applied to these HDF- 
N flux density results because in many cases the spectral 
shapes and redshifts of objects are not known adequately 
enough to apply any correctio n with confidence. Ad dition- 
ally previous studies, such as bv lAppleton et al. (2004), have 
shown that the application of k-corrections to individual 
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Figure 9. The luminosities at 1.4 GHz and 24 fim for 224 24 /im Spitzer sources with redshifts and significant radio flux density from 
the HDF-N (black triangles) and 50 non-AGN 2A/im sources with redshifts catalogued by Norris et al. (2006) (green triangles). The 
luminosities shown have not been k-corrected. All 1.4 GHz luminosities have measured within an aperture of radius 1'.'5 centred upon 
the Spitzer MIPS 24 (im source position. Note the few cases where HDF-N sources with apparently low radio luminosity with respect to 
MIR luminosity are large nearby spirals with an optical extent larger than the 3". 



sources below z=l only introduces a small change to the 
average value of q24 as a function of redshift. Of the sources 
with known redshifts in our sample 78 percent have redshifts 
less than 1, and thus will be minimally affected. Within 
Figs.[5]and[7]the values of q24 for individual sources are plot- 
ted against their flux densities at each wavelength. These 
sources sample a wide range of luminosities and redshifts 
and hence the systematic effects resulting from not apply- 
ing a fc-correction will effect individual sources but will cause 
only minimal changes in any overall trend of the sample as 
a function of flux density. However, the effects of the ap- 
plication of a ^-correction to these data do warrant further 
discussion. 

Assuming that the 24 /im sources identified in this sam- 
ple are predominately powered by star-formation, at an 
observed radio frequency of 1.4 GHz the emission will be 
dominated by non-thermal synchr otron emission from su- 
pernovae and supernova remnants l|Muxlow et al.ll2007h . As 
such the non-thermal radio signal will be boosted by a factor 
of (1+z) ' 7 , assuming a synchrotron power law (Soc u~ ' 7 ). 
This ^-correction will become incorrect for the highest red- 
shift sources as the rest-frame radiation will be shifted 
to higher radio frequencies where the emission from other 
mechanisms such as thermal free-free emission will become 
increasingly important. For example a source at redshift of 
i52.5 the radio emission observed at 1.4 GHz will have been 
emitted at a rest-frame frequency of si 5 GHz at which point 
an increasing proportion of the radio emission from a star- 
forming galaxy will arise from thermal processes resulting in 
a reduction in the radio spectral slope and hence ^-correction 
that should be required. 

Accurately fc-correcting the Spitzer 24 (im flux densities 
requires detailed knowledge of the mid-IR spectral template 



of the individual sources. In particular, for sources at red- 
shifts > 1 significant spectral features from P AHs and sili- 
cates observed in local starburst galaxies (e.g. iBrandl et all 
2006) are shifted into the Spitzer 24 fxm band which can re- 
sult in large additional fc-correction factors. The inclusion of 
these line features can result in an increase of a factor of 2 
or more in S24/im compared with commonly used sta rburst 
spectral templates l|Fadda et aZj2006l : lYan et qZ J2007f ). This 
is equivalent to a change in q24 of approximately 0.3. The 
general effect of the MIR fc-correction, resulting from a spec- 
tral index of a starburst galaxy excluding line emission, will 
be to boost S24 M m as a function of redshift by an amount 
that will depend upon spectral slope of the galaxy template 
used. For a luminous starburst this will result in a boost - 
ing with a MIR spectral slope a\ 5 ~ 2 (|Yan et qZ.ll2007f ). 
As such the effect of this k-correction applied to the 24 /im 
data will be larger than that applied to the 1.4 GHz data, 
hence the overall effect will be to increase the value of q24 
as a function of redshift. 

Uncorrected q24 is plotted against redshift for the 
sources with known spectroscopic or photometric redshifts 
in Fig[5] Applying ^-corrections to these data points will re- 
sult in a small increase in q24 for the highest redshift sources 
plotted. Within the scatter of the small sample of sources 
plotted in Fig[S]no significant trend is observed as a function 
of re dshift to z~3. This is consistent wi th previous studies 
(e.g. lGarrettll200l ; lAppleton etai\\2004 ) ■ 



4-2.1 Potential contamination from AGN 

In addition to instrumental and analysis effects outlined 
above it is possible that the flux densities of sources recorded 
in this study may include some proportion of contamina- 
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tion by embedded AGN l|Richards et aZj|2007f ). The effect 
of the inclusion of sources with significant AGN emission is 
demonstrated by t he AGN sources from the SWIRE survey 
|Norris et al\\200(h which are overlaid as stars in Figs. [2] [5] 
and [Jj In general, sources with a significant AGN contri- 
bution show an excess of radio emission compared to their 
infrared luminosity, hence will deviate from the IR-radio cor- 
relation. Consequently the result of any AGN contamination 
within the HDF-N Spitzer sample analysed here will be to 
reduce q24. However, the flux densities of the HDF-N sam- 
ple (see Fig. [2} show only a very small number of sources 
with excessive radio flux density which deviates significantly 
from the IR-radio correlation extrapolated from high flux 
density samples. This is to be expected since it is now gen- 
erally accepted that below 100 /iJy the vast majority of ra- 
dio sources are predominantly powered by star-formation 
llRichardsll200d;lGruppioni et aZ.ll2003l ; ISevmour et afj|2004t 
iMuxlow et a/.l 12005). Additionally this sample of sources 
is selected at 24 fim rather than at radio wavelengths and 
hence will preferentially select star-forming sources rather 
than radio AGN. Thus, throughout this paper no direct dis- 
tinction has been made between possible AGN and non- 
AGN star-forming galaxies in these faint Spitzer selected 
HDF-N radio sources. However, a small AGN contamina- 
tion within the sources sampled cannot be eliminated. The 
consequences of any AGN contamination may, to some de- 
gree, offer an explana tion for the lower v alues of q24 derived 
here (see for example iDaddi et aZ.ll2007h . 



4.3 Causes and implications of an evolution in the 
IR-radio correlation 

In the previous section a small trend in q24 as a function 
of 24 /im flux density has been tentatively identified. This 
trend is small and only visible in statistically averaged deep 
radio observations. The much deeper radio surveys that will 
be made using the newly enhanced radio instruments such 
as e-MERLIN and the EVLA will constrain the faint end of 
the IR-radio correlation, confirming or refuting any possible 
deviations. 

Whilst there are several observational and analysis ef- 
fects which may contribute to the observed trend in the IR- 
radio correlation (outlined in the previous sections) there 
are also physical mechanisms which may cause an effect 
similar to that tentatively observed. Both the IR and ra- 
dio emission from star-forming galaxies are indirectly linked 
with the rate of ongoing star-formation; simplistically the 
non-thermal radio emission originates in accelerated charged 
particles in supernovae and supernova remnants and the IR 
emission results from the reprocessing by dust of light from 
hot young stars. This inferred connection is generally con- 
sidered to result in the IR-correlation. However since both 
the radio and IR emission are linked with star-formation via 
complex physical processes it is remarkable that the lumi- 
nosity in these two bands correlates so accurate ly ov e r such 
a large range of luminosities. A recent study by iBelll (|2003T ) 
has provided a detailed investigation of the relationship of 
the radio and total infrared emission of galaxies and their 
star-formation rates. In this work several competing physi- 
cal factors which will affect the IR and radio emission from 
star-forming galaxies are considered which at low luminosi- 



ties potentially cause a deviation in the IR-radio correlation 
similar to that which is very tentatively reported here. 

If we assume that the non-thermal radio emission 
from star-forming galaxies is directly proportional to the 
star- formation rate (i.e. radio oc SFR), then the system- 
atic reduction in q24 at low 24 fim flux density implies 
that the 24 /xm emissi on is no longer efficiently tracing 
the SFR of the galaxy fe cll 2003) . Observational ev i dence 



(e.g. IWang fc Heckmanl Il99rj; lAdelberger fc Steidell l200d: 



Hopk ins et al l200ll ; iBuat et all2002h has shown that galaxy 



luminosity and dust opacity to UV and Ha emission are cor- 
related, demonstrating that lower luminosity galaxies have 
substantially less dust absorption and reddening than high 
luminosity sources. These studies have also shown that the 
ratio of far-UV to IR light is significantly lower in less lu- 
minous sources, implying that whilst for luminous galaxies 
the IR provides a good tracer of star-formation at low lumi- 
nosities much of the far-UV emissi on from star-formation is 
not reprocessed as IR emiss ion (see IWang fc Heckmanlll996l ; 
IBuat et aZ.ll2002i ; lBell|[2003l . for further discussion). The con- 
sequences of the reduced dust opacity in lower luminosity, 
presumable smaller, sources will be a reduction in the effi- 
ciency of reprocessing of starlight into IR emission and hence 
will result in an observable decrease in the ratio of IR-to- 
radio emission from sources at low luminosities. Of course 
it should be reiterated that this makes the assumption that 
the radio emission is a 'perfect' tracer of the star-formation 
which may equally not be true. 

It has been suggested that the synchrotron radio emis- 
sion from low- luminosity galaxies can also be suppressed, 
implying that it is not a true trace r of the SFR (|Klein et all 
1 19841 . Il99ll ; iPrice fc Duricl ll992T l. Whilst the connection 
between non-thermal radio emission and star-formation 
is complex, it has been discussed, both on the basis of 
models of the physical mechanisms and obse rvational ev- 
idence, by many auth o rs (see for example IVolk I 1 19891 ; 
IChi fc Wolfendaleill990l ; iHelou fc Bicavl Il993h . A physical 
explanation for the suppression of the non-thermal radio 
emission in low-luminosity galaxies is that the accelerating 
electrons can escape these sources because they have less ef- 
ficient c osmic-ray confinement th an in larger, more luminous 
sources (|Chi fc Wolfendalelll990r ). 

Radio surveys of IR sources prior to the very deep re- 
sults presented here have not detected any significant de- 
viations in the lin e arity of the radio-IR correlation (e.g . 
Condon et o71ll99ll : lYun et aZ-lkoOll ; lAppleton eTal]\2004 ). 



Belli (|2003T ) argue that this observed linearity at relatively 
low luminosities implies that both the IR and non-thermal 
radio emission do not adequately trace the star-formation 
in low luminosity sources. As such the observed linearity of 
the IR-radio correlation arises because the two independent 
physical mechanisms by which the emission in these bands 
originates both result in an underestimate of the SFR at 
low luminosities by a similar amount, thus preserving the 
observed linear na ture of the correlation in a conspirato- 
rial way l|BeHll2003T ). However, these two competing physical 
trends whilst similar in magnitude will diverge for very low 
luminosity sources resulting in a reduction in the IR-to- radio 
ratio . This effect can be seen in (figures 3 and 8 from IBelll 
2003) and is tentatively hinted at in the observation pre- 
sented here (Fig. [5}. 

Not withstanding the factors, discussed in the previous 
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section which may effect this tentative detection of a small 
deviation in the gradient of the IR-radio correlation, these 
obser vational res ults are broadly consistent with the mod- 
els of Bell i|2003h . If this observed evolution of the IR-radio 
correlation at low flux densities is real it potentially has 
significant implications on the use of either IR or radio lu- 
minosities for the derivation of star formation rates of faint 
galaxies in future deep surveys. 

4.4 Future observational tests 

Future observations will be required to confirm this tentative 
deviation of the IR-to-radio correlation at these extremely 
low flux densities. These results may be tested further by i) 
deeper radio and infrared observations which detect many 
faint individual sources at high significance, ii) significantly 
larger area deep surveys that allow robust statistical detec- 
tions of faint sources to be made, or Hi) the separation and 
analysis of the various competing factors which contribute 
to this effect. 

The stacked radio images of extremely faint IR sources 
presented in Fig. [3] provide a tantalising glimpse of the im- 
ages which will be routinely available in the near future. The 
next-generation of radio interferometers, such as e-ME RLIN 
(|Garrington et qZ. .11200-3 ) and the EVLA jPerlev Il200ch . will 
provide images many tens of times more sensitive than can 
currently be made. Future deep radio observations using e- 
MERLIN and the EVLA with integration times compara- 
ble to those presented here will produce images with sub- 
microjy noise levels. Surveys using these instruments will 
easily be able to detect all of the IR sources so far found in 
the deepest Spitzer observations with very high confidence, 
and the high sensitivity of these instruments will allow much 
larger areas of micro Jy radio sky to be surveyed. Such high- 
sensitivity and wide-field observations will constrain any low 
flux density deviations in the radio-IR correlation. 

The radio observations presented here primarily trace 
the non-thermal synchrotron radio emission of galaxies with 
flux densities of a few tens of /iJy. At 1.4 GHz ~90 per- 
cent of the emission of a typical ~L* spiral galaxy is non- 
thermal emission with the remainder com ing from thermal 
bremsstrahlung processes (|Condonlll992T ). Whereas the ra- 
dio synchrotron, as stated previously is indirectly related 
to star-formation, the bremsstrahlung emission is a more 
direct tracer of the ongoing star-formation since it arises 
from thermal radio emission from gas ionised by hot, young 
stars. At higher radio frequencies (>10GHz) thermal ra- 
dio e mission will begin to dominate over s ynchrotron emis- 
sion l|Condon |[l992l ; iPrice fc Duric 1 1 19921 ) tending to flat- 
ten the radio spectrum of star-forming galaxies. Determin- 
ing the relative fraction of thermal and non-thermal radio 
emission in different sources and at different frequencies is 
non-trivial requiring multi-frequency measurements of the 
radio flux densities in order to determine the slope of non- 
thermal emission. Such observations will be more feasible 
using future extremely deep radio observations planned with 
the new generation of wide-bandwidth radio interferometers, 
such as e-MERLIN and the EVLA, which will allow sensitive 
mult i- frequency radio observations of the distant Universe to 
be made within reasonable observation times. Using obser- 
vations made with these enhanced centimetre instruments 
along with measurements made using longer wavelength ra- 



dio interferometers like the Giant Metrewave Radio Tele- 
scope (GMRT) and LOw Frequency ARray (LOFAR) it will 
be possible to constrain the radio spectral energy distribu- 
tion of many star-forming galaxies hence allowing the ther- 
mal and non-thermal radio emission to be separated. By 
disentangling the thermal and non-thermal radio luminosi- 
ties of many distant star-forming galaxies it will be possible 
to verify the expected suppression of IR emission from low 
luminosity sources relative to the thermal radio emission 
which is a more direct, extinction free, tracer of the star- 
formation. 



5 CONCLUSIONS 

Using one of the deepest high- resolution 1.4 GHz observa- 
tions made to date, in conjunction with deep 24 /jm Spitzer 
source catalogues from GOODS, we have investigated the 
microjy radio counterparts of faint MIR sources. These ob- 
servations confirm that the microjy radio source population 
follow the MIR-radio correlation and extend this correlation 
by several orders of magnitude to very low flux densities and 
luminosities, and out to moderate redshifts. This extension 
of the MIR-radio correlation confirms that the majority of 
these extremely faint radio and 24 /jm sources are predomi- 
nantly powered by star-formation with little AGN contami- 
nation. 

Statistically stacking the radio emission from many tens 
of faint 24 fim sources has been used to characterise the size 
and nature of the radio emission from very faint IR galaxies 
well below the nominal radio sensitivity of these data. Using 
these methods the MIR-radio correlation has been further 
extended and a tentative deviation in this correlation at very 
low 24 pirn flux densities has been identified. 
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